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To resolve the puzzle of the proton size raised from the recent result of muonic 
hydrogen Lamb shift, De Rujula has proposed that a large value of the third Zemach 
■ moment (r~)(2) of the proton to be the solution. His suggestion has been criticized 

(N : 

q^. by many groups based on the ep scattering data at low Q 2 regime. However, if there 

\ is a "thorn" or "lump" in the electric form factor of the proton Ge(Q 2 ) at extremely 

low Q 2 regime, then the third Zemach moment (fp)(2) would be as large as De Rujula 
suggested. In this article, we show that the existence of such a "thorn" or "lump" 
^ , has not been completely excluded, although tightly restricted, by the current data of 

J|h ; ep elastic scattering. We also suggest a more sophisticated global fitting procedure 

of Ge(Q 2 ) for the future fitting. 
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The issue of the charge radius of the proton has attracted a lot of attention, since the 
charge radius extracted from the Lamb shift of muonic hydrogen has been reported to be 
0. 84184(67) fm [1]. This result is significantly smaller than the previous value of CODATA 
^ . |2| , J (r^) (CODATA) =0. 8768(69) fm and the one extracted from ep elastic scattering data 

a ' i — 



I. INTRODUCTION 



(r 2 )(ep)=0. 879(5) s tat(4) sys t(2) mode i (A) group fm J3j. De Rujula [4J, |5] has pointed out that 
the small value of the charge radius reported in |l| based on the assumption that the electric 
from factor of the proton Ge(Q 2 ) is the dipole form. Hence the original QED formula, 

L th (meV) = 209.9779 - 5.2262(rJ) + 0.00913(r^) (2) , (1) 

is reduced into 

L th (meV) = 209.9779 - 5.2262(rJ) + 0.0347(rJ) 3/2 , (2) 

because 

m m ?~[(rl)]\ (3) 
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when Ge(Q 2 ) is the dipole form. (Note that in the above equations, the units of (r 2 ) and 
(fp)(2) are fm 2 and fm 3 , respectively.) Since the experimental result is L e:rp =206.2949 ± 
0.0032 meV, accordingly they concluded that the value of the charge radius is 0.84184 fm 
However, De Rujula has argued that there is no reason to believe Ge{Q 2 ) to be the 
dipole form. Instead he suggested that the proton may own a large third Zemach moment 
about 36.59 fm 3 , which is about fifteen times larger than the value from Eq.(|3]). If so the 
value of the charge radius extracted from Eq.([T]) will agree with the CODATA value well. 
Furthermore he has employed some toy model of the GeIQ 2 ) to obtain (r 3 )( 2 )=36.59 fm 3 . 
By this way one is able to resolve the proton size puzzle |4j, |5 ] . 

The other attempts to resolve this puzzle, for example, to recalculate the polarizability 
contribution [6] or to estimate the non-perturbative effect [7], and to test the possibility 
of the existence of the new particle between the proton and the muon have been not 
very successful so far. The new corrections they have found are usually too small (The only 
exception is so-called off-mass-shell effect advocated by [9]). Therefore the simple solution 
suggested by De Rujula seems to be worthy of further investigation. 

However, the proposal of De Rujula has been severely criticized by several groups 



12]. First, the toy model used by De Rujula has been indeed ruled out by the recent 



experimental data. Furthermore, those groups have argued that such a large value of the 
third Zemach moment cannot accommodate the current data of ep elastic scattering. What 
they did, instead, is to adopt several widely used parametrizations of Ge{Q 2 ) to calculate 
the correspondent third Zemach moment (0(2) and presented them to be far smaller values 
than the one obtained by De Rujula. At first glance, this objection looks very convincing. 
However, as De Rujula already pointed out, the value of (r 3 ,)^) is extremely sensitive to the 
behaviour of Ge{Q 2 ) in very low Q 2 regime. Because there is no data between Q 2 = to 
Q 2 = Q 2 min - Therefore he argued the slim possibility of large third Zemach moment may 
not completely excluded yet [5]. 

But one can provide a counterargument as follows: the extrapolation of Ge{Q 2 ) between 
Q 2 =0 and Q 2 =Q 2 nin should be very reliable. Because the value of Ge at Q 2 = must be 
one due to the fact that the electric charge of the proton is +e, and its derivative - at 

Q 2 =0 is also severely constrained by the CODATA value of (r 2 ). Thus the extrapolation of 
Ge{Q 2 ) from Q 2 =Q 2 nin to Q 2 =0 is supposed to be adequate enough to determine the value 
of (r 3 )(2). However this counterargument has one loophole. If there appears a "thorn" or 
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' dip " in Ge{Q 2 ) between Q 2 =0 and Q 2 =Q min then the parametrizations previously used 
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14J will no longer be able to produce an accurate value of (?"p)(2)- Naturally one should 
ask that whether there exists a Ge{Q 2 ) with "thorn" or "lump" which can generate a large 
( r p)(2)> an d at the same time, accommodate the existing ep scattering data. In particular, 
recently a measurement of the cross section of the elastic ep scattering has been carried out 
at Mainz, rang ed from Q 2 = 0.004 GeV 2 to 1 GeV 2 with the statistical errors below 0.2% 
jsj. Their data has shown no sign of any "bump". It makes any attempt to obtain large 
(r 3 )( 2 ) by adding "thorn" at Ge(Q 2 ) to be very difficult. But in this article we will explicitly 
show that such a task is indeed difficult but not totally impossible. 

The outline of this article is as follows. We first review the relationship between the 
third Zemach moment and the electric form factor Ge{Q 2 )- Next we combine the "thorn" 
and some parametrizations of Ge{Q 2 ) to calculate the third Zemach moment and learn the 
relation between the height, width and peak position of the "thorn" and the third Zemach 
moment of the proton. Then we explicitly show that one can combine our ansatz of "thorn" 
with the inverse-polynomial fit used in |3[ to obtain a large third Zemach moment as De 
Riijula has suggested. At the same time the combined ansatz deviates from the original 
inverse-polynomial fit less than 0.2%. Finally we present our conclusions and outlooks. 



II. ZEMACH MOMENT AND THE ELECTRIC FORM FACTORS 

The conventional proton charge density is defined as the Fourier transform of the electric 
form factor Ge{Q 2 ) in the Breit frame, 

P P (r) = J (^^^(q)- (4) 

Here Q 2 = —q 2 = \q \ 2 — g 2 which is equal to \q \ 2 in the Breit frame. We use the notation 
q = \q\. The following quantities are defined as 

(r»> = J d 3 rr n Pp (r). (5) 

From this definition one can easily deduce that Ge(0) = 1 because the 0-th moment (r°)=l 
and rfG j^2* - [q 2 =o = ~\( r p)- On the other hand the n-th. Zemach moment is defined as 

(r n )(2) = J d 3 rr n p 2 (r), (6) 
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where P2( r ) is defined as 



P 2 (r) = J d 3 r'p p (r')p p (r'-r) = J ^ e -**G|(q) 



(7) 



After some algebra one obtains the following result [l5|, 



(2) 
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q 2 , , 



^(q)-^(r )-l 



(8) 



It is obvious that the third Zemach moment of the proton is dominated by the Ge{Q 2 ) at 
very low Q 2 . The crucial issue here is whether there exists one form of Ge{Q 2 ) which is able 
to produce large (t%)(2) and at the same time accommodate the current data of ep elastic 
scattering. 



III. RELATION BETWEEN THE THORN IN Ge AND THE THIRD ZEMACH 

MOMENT 

In this section we assume that there is some "thorn" or "lump" appearing in the Ge{Q 2 ) 
in the very low Q 2 regime. We expect such a pathological structure to produce a large third 
Zemach moment. Here we express the electric form factor as follows, 

Ge(Q 2 ) = G { *\Q 2 ) + AGe(Q 2 ), (9) 

where G£\Q 2 ) is some parametrization from the global fitting of the ep scattering data. On 
the other hand AGe{Q 2 ) denotes the "thorn" on the electric form factor. Naively, one may 
think that it is easier to simply add a triangle function with the height H and the width W, 
whose peak is located at Q peak - However such a choice will cause a serious problem. One can 
calculate the associated charge density Ap(r) by making the Fourier transform of AGe{Q 2 ), 
then calculating its contribution to (r 2 ). However, if the triangle function is chosen then its 
corresponding A(r 2 ) calculated by 

A(r 2 p ) = J d 3 rAp(r)r 2 , (10) 

is actually divergent! It is due to the fact of the corresponding Ap(r) actually converges 
slower than 1/r 4 . Hence one has to make judicious choice of the "thorn" function so that 
(r 2 ) can be kept finite. On the other hand, here we still want to employ the widely used 
parametrizations of Ge(Q 2 ) whose value at the Q 2 = have been fixed. As a result 
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AGe(Q 2 = 0) and d ^ 2 E (Q 2 = 0) both have to be negligible. Moreover the influence of 
AGe(Q 2 ) has to able to be ignored when Q 2 > Q 2 min - One needs figure out some function 
form satisfying the above criteria. Here we present our choice as follows, 



AG E {Q 2 ) = K\ exp 



(Q 2 - K 2 ) 2 



(11) 



Here K\ is dimensionless and the unit for K 2 and K% is GeV. Ki, K 2 and K3 denote the 
height, the position of the peak and the width, respectively. 

To explain the result of muonic hydrogen Lamb shift one needs show that the following 
quantity 

AL(meV) = L theory - L exp = 209.9779 - 5.2262(r 2 ) + 0.00913(r 3 ) {2) - 206.2949, (12) 

to be smaller than the experimental uncertainty 3 x 10~ 3 meV. If we choose the parametriza- 
tions of [kJ or [l^] as our G^\ it is easy to pick up several parameter sets of -^1,2,3 to satisfy 
all criteria. We list our parameter sets and their corresponding values of (r 2 ) and (r 3 )( 2 ) 
in Table (I). One may wonder the values of (r 2 ) is somehow too small compared with the 
CODATA value: (r 2 )(CODATA)=0.753 fm 2 . The reason for it is because the parametriza- 
tions we used are the results of global fitting and their value of (r 2 ) are somehow small. For 
example, (r 2 )(Albrico)=0.750 fm 2 and (r 2 ) (Kelly) =0.744 fm 2 . 





r (R) 

E 


Ki 


K 2 (GeV) 


K 3 (GeV) 


AL(meV) 


(r 2 ) (fm 2 ) 


(r 3 ) (2) (fm 3 ) 


I 


Alberico 


0.119185 


0.08 


0.0447214 


1.7xl0" 4 


0.745137 


23.138 


II 


Alberico 


0.0139929 


0.08 


0.053183 


1.5xl0" 4 


0.717153 


7.11973 


III 


Alberico 


0.283648 


0.10 


0.053183 


-2.85xl0" 4 


0.747693 


24.5962 


IV 


Alberico 


0.139056 


0.10 


0.0588566 


8.56xl0~ 6 


0.735336 


17.5272 


V 


Alberico 


0.720091 


0.12 


0.053183 


1.79xl0~ 4 


0.748308 


24.9536 


VI 


Kelly 


0.130982 


0.08 


0.0422949 


1.06xl0~ 6 


0.742016 


21.3496 


VII 


Kelly 


0.101611 


0.08 


0.0447214 


-9.31xl0" 6 


0.739647 


19.9926 


VIII 


Kelly 


0.243554 


0.10 


0.053183 


1.06xl0" 4 


0.741824 


21.2389 


IX 


Kelly 


0.118405 


0.10 


0.0588566 


-8.5xl0~ 4 


0.731308 


15.2113 


X 


Kelly 


0.627337 


0.12 


0.053183 


4.6xl0~ 6 


0.742354 


21.5437 



TABLE I: Our chosen parameter sets and the values of their corresponding (r 2 ) and (rp)r 2 y 
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Unfortunately the above results have all been excluded by the recent Mainz low Q 2 
dataQ. Nevertheless, we have observed several important facts from the Fig (1). First, if 
we make the position of peak, K 2 , more close to the Q 2 =0, the height of the peak will be 
smaller with the same width. However, if K 2 becomes too small, it will produce relatively 
large A(r 2 ), which is negative, thus the resultant (r 2 ) is much smaller than the CODATA 
value. The second important fact is as follows. With the same K 2 , the height K\ decreases 
as the width K 3 increases, i.e., when the peak is less sharp and the height becomes smaller. 
Thirdly, we also find that the result is not very sensitive to the choice of the Gg(Q 2 ) 
parametrizations as shown by the Table (I). These facts will instruct us to construct more 
realistic ansatz of Ge{Q 2 ) as shown in the next section. 

GE 

1.8 r 

1.6 - 
1.4 - 
1.2 - 
1.0'- 
0.8 - 



FIG. 1: The unit in x axes is GeV 2 . (Left)The G E {Q 2 ) in (II)(solid line), (III) (dashed line) and 
(V) (dotted line). Three curves correspond to the same i^3=0.053 GeV but their values of K 2 are 
0.08 GeV, 0.10 GeV and 0.12 GeV respectively. (Middle)The G E (Q 2 ) in (I) (dashed line), (II) (solid 
line) with the common value of K 2 =0.08 GeV. But their corresponding values of K3 are 0.0447 
GeV and 0.0532 GeV, respectively, (right) The G E (Q 2 ) in (III) (dashed line), (IV) (solid line) with 
the common value of i^2=0.10 GeV. But their corresponding values of K3 are 0.0532 GeV and 
0.0589 GeV, respectively. 




IV. INVERSE-POLYNOMIAL FIT WITH A "THORN" OF G E 



We have learned how to increase the third Zemach moment by adding the "thorn" in 
the previous section. Here we need construct a parametrization to accommodate the recent 
Mainz low Q 2 data with Q 2 min — 0.004 GeV 2 with the uncertainty below 0.2%. In this section 
we show that one is able to combine our ansatz of "thorn" with the inverse-polynomial fit 



used 



in (3 



161 ] to make the difference between L th and L exp to be smaller than 3 x 10 3 meV. 
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At the same time the combined ansatz deviates from the original inverse-polynomial fit less 
than 0.2%. The inverse-polynomial fit has been used in {2I [if]]. Its explicit form is given as 

1 



Qinv—poly ^q2 



1 + Ei=i a iQ 



2i 



(13) 



with 



ai = 3.3615, a 2 = -3.0343, a 3 = 29.6677, a 4 = -85.6169, 
a 5 = 130.7053, a 6 = -101.5145, a 7 = 34.2926. 

yinv— poly / 



(14) 



Note that G l ^ v ~ pmy (Q 2 ) generates (r 3 )( 2 )=2. 96667 fm 3 . To combine this fit with the ansatz 
in Eq. ljTTj) . one needs to guarantee that G E (0) = 1. Hence we modify Q™ v -v ol y j n ^ the 
following one, 

GT d (Q 2 ) = (1 - AG E (0))G™- poly (Q 2 ) + AG E (Q 2 ). (15) 

Employing the ansatz in Eq. ( 1T51) , one can make AL to be smaller than 3 x 10~ 3 meV with the 
chosen parameters listed in the Table (II). The value of {r 2 ) is a little larger than CODATA 
value but still be reasonable. To accommodate the very precise Mainz low Q 2 data, it is 
necessary to make K 2 very small. The width K 3 is about only half of the value used in the 
previous section. The height K\ is only few percents of the ones in Table (I). However such a 
"lump" generates a very large (r 3 )( 2 ). It shows that the third Zemach moment is extremely 
sensitive to the detail of the electric form factor Ge{Q 2 ) at very low Q 2 regime. 





K 2 (GeV) 


K 3 (GeV) 


AL(meV) 


{r 2 ) (fm 2 ) 


(r 3 ) (2) (fm 3 ) 


-0.0016962 


0.001 


0.0221336 


5.95xl0- 4 


0.787349 


47.3646 



TABLE II: Our chosen parameters and the values of (r 2 ) and (r 3 )( 2 ). 



Moreover we define the following quantity to characterize the difference of our modified 
fit and the original inverse-polynomial fit. 



R(Q 2 



G™ d {Q 2 ) - G^ v - poly {Q 2 ) 



y^iinv—poly 



(Q 2 



(16) 



The left panel of Fig. (2) shows that the value of R{Q 2 ) is ranged from to 0.17%. It 
is due to the fact AG E {Q 2 = 0) = —1.6 x 10~ 3 . When Q 2 is large enough, AG E ~ and 



G mod {Q 2 ) ~ (1 - AG E (0))G 



inv—poly 



(Q 2 ). One can also observe the curve of G mod (Q 2 ). It is 
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FIG. 2: (left) The value of R(Q 2 ) defined in the text, (right) The curve of G^ od (Q 2 ). 

a very smooth "lump" hidden at the extreme low Q 2 regime. The shape of this "lump" is 
depicted in the right panel of Fig. (2). It is very smooth as one can observe from the plot. 

One may frown on our result here and argue that our result cannot accommodate the 
Mainz data satisfactorily. Indeed the ansatz in Eq. (fT5|) is quite a simple way to add a "lump" 
to the existing fit. However, we emphasize that even with such a simple ansatz one can still 
embed a "lump" at Ge(Q 2 ) and produce a large (r^)^). We believe it to be promising to 
improve the result with better agreement with the data by employing more sophisticated 



ansatz instead of the one we used here. We leave it for our future publication [17 ]. 



V. CONCLUSION AND OUTLOOK 



In this article, we show that the third Zemach moment becomes large if there is a "thorn" 
or "lump" at very low Q 2 regime. Furthermore, our study show that it is possible to construct 
a parametrization of Ge{Q 2 ) which can accommodate the existent ep elastic scattering data 
and, at the same time, generate a large (fp)(2) to explain the Lamb shift of the muonic 
hydrogen. 

In this work we limit ourselves to combine the existent parametrizations of G l E lv ~ poly (Q 2 ) 
and a simple ansatz denoting the " thorn" by a very simple way in Eq. ( TT5T) . In principle one 
should use the following ansatz, 

ai + a 2 Q 2 



G E (Q< 



+ bi exp 



-(Q 2 -b 2 ) 2 



1 + a 3 Q 2 + a 4 Q 4 + a 6 Q 6 
to fit the ep scattering data globally. There are two relations between those parameters, 



(17) 



1 = a\ + b\ exp 



-6f- 



(18) 
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(r 2 p )(CODATA) 
6 



— a 2 + a 3 + 



1h x b\ 



(19) 



Here we have four new parameters a\ and &i_3 with several constrains such as Eq. ( fl8|) 
and Eq.f ll9p . The third constraint is the resultant result of Eq.([H|) has to be around 36 fm 3 . 
Using the above parametrization, it is likely that one can pick up a suitable parameter set 
to accommodate the existent ep elastic scattering data. The result by default can explain 
the Lamb shift of both electronic and muonic hydrogen. We leave this task for our future 
publication [171 ]. 

Although phenomenologically a large third Zemach moment is possible, nevertheless, 
there are still many challenges from theory side. For example, the very low Q 2 behaviour 
of Ge{Q 2 ) is supposed to be dominated by the chiral physics. Namely the pion cloud plays 
the crucial roles in the low energy regime and one can apply Chiral Perturbation Theory (x 
PT) to calculate the electric form factor there [18|. The x PT result of (?"p)(2) is about 2 — 3 
in 3 , which is much smaller than ours. We also notice the most recent estimate made by 



19j , their conclusion disagrees with us. The reason is they insist to adopt the smooth p(r) 
which guarantees (r n ) to be always finite. We only require the convergence of Eq.([H]) and 
Eq.f JTQ|) only. These issues all remain open and need further studies. 
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